arXiv: 1501.02786vl [cond-mat.mes-hall] 12 Jan 2015 


First-order correlation function of the stream of single-electron wave-packets 

Michael Moskalets* 

Department of Metal and Semiconductor Physics, 

NTU “Kharkiv Polytechnic Institute”, 61002 Kharkiv, Ukraine 
(Dated: January 13, 2015) 

The first-order correlation function, which is accessible experimentally, contains all essential infor¬ 
mation about the state of the system of non-interacting electrons. Here I discuss how this function 
can be used to answer the question whether the state of a periodic stream of single-electron wave- 
packets is a multi-particle state or it is the product of single-particle states. In the latter case the 
correlation function is expected to be factorizable while in the former case it is not. As an example 
I consider a train of Lorentzian in shape single-electron excitations, levitons. I demonstrate that the 
correlation function in time domain is factorizable or not depending on whether the wave-packets 
are separated or overlapping. In contrast, the correlation function in energy domain is always 
factorizable and thus cannot be used to distinguish single- and multi-particle states. 
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I. INTRODUCTION 

Recently Jullien et al} have reported a measurement of 
a wave function of a single electron traveling in a ballistic 
conductor. This groundbreaking experiment opens the 
route to the control and characterization of single- and 
few- electron states injected on demand into solid state 
circuits. 

The experiment of Ref. 1 was done with a train of elec¬ 
trons, each prepared in a well-controlled single-electron 
state called a leviton. Such a train is also appropriate 
to address another problem important for a few-electron 
state engineering. Namely, it is important to demon¬ 
strate directly that a multi-particle state is not formed 
in the stream of levitons. Such a multi-particle state 
should be necessarily formed if the distance between the 
levitons would become comparable with the width of a 
leviton, see Fig. 1. So, merely by varying the rate at 
which the levitons are generated, one can tune the state 
of a stream from a multi-particle state to (the product 
of) single-particle states. 

In quantum optics^^ the single- and multi-photon 
states are distinguished with the help of the intensity- 
intensity correlation measurement^^d4. photon flux 
is divided by a beam splitter and the resulting fluxes are 
directed to two single-photon detectors. The absence of 
coincident detections proves that a multi-photon state is 
not present in the flux. This technique is customarily 
used to characterize single-photon sources.^® 

In quantum coherent electronics the analogue of an 
optic beam splitter is a quantum point contact, 
which divides an electric current into two beams, the re¬ 
flected and the transmitted one. Thus the current cross¬ 
correlation measurement is expected to be an analogue 
of a coincidence detection measurement in optics. How¬ 
ever this is not the case. The reason is that there are 
no efficient single-electron detectors available at the mo¬ 
ment. Therefore, what is measured in electronics is a 
current noise, the product of fluctuations of reflected and 
transmitted currents averaged over a long time. Apart 
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FIG. 1. (Color online) The periodic stream of single-electron 
wave-packets: If the width of an individual wave-packet is 
smaller than the period, i.e., the distance between subsequent 
wave-packets, then all electrons are distinguishable by the pe¬ 
riod they occupy and no a multi-particle state is formed. In 
particular during each period only one electron can be de¬ 
tected. If in contrast the width of a wave-packet is comparable 
with the period then more then one electron can be detected 
during a period, though the mean number of detected parti¬ 
cles remains the same. Due to the overlap between different 
wave-packets the multi-particle state is formed. 


of the thermal fluctuations vanishing at low tempera¬ 
tures, the noise contains a part, known as the shot noise, 
which relies on the quantization of charge and exists due 
to the fact that an electron can be either transmitted 
or reflected but not both.^® In the absence of electron- 
electron correlations the shot noise is proportional to the 
mean number of particles in a stream per unit time,^®^^® 
since each particle is scattered independently. This fact 
was used in Refs. 8 and 24 to count directly the number 
of particles emitted in each cycle. Another approach to 
characterize a single-electron emission regime is to mea¬ 
sure a finite frequency noise, which also provides infor¬ 
mation on a short-time dynamics of the source. 

If the distance between particles within a periodic 
stream would decrease such that different particles start 
to overlap, see Fig. 1, the mean number of particles 
per period would not change while a multi-particle state 
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would emerge. This results in fluctuations of the num¬ 
ber of particles measured during a given period. The 
shot noise, however, cannot distinguish single- and multi¬ 
particle states as long as the mean number of particles 
per period is the same. 

The aim of this work is to demonstrate that another 
measurable, namely the first-order correlation function, 
can be used to make the necessary distinction. In Refs. 29 
and 30 it was shown that the correlation function in 
time domain is directly accessible via a time resolved 
interference current, which is measured at the exit of 
an electronic Mach-Zehnder interferometer^^. Another 
approach to access a correlation function but in energy 
domain is a single-electron quantum tomography.^^^^® It 
relies on a zero-frequency noise measurement in the Han- 
bury Brown and Twiss set-up^®^^®, where the source of 
interest is in one input channel while the other input 
channel is fed by a weak ac voltage. Such an electronic 
quantum tomography was realized in Ref. 1. 

The paper is organized as follows: In Sec. II some ba¬ 
sic properties of fermionic correlation functions are pre¬ 
sented and the energy representation in the case of a 
periodic flux is introduced. In Sec. Ill the correlation 
function of a periodic train of electron-like excitations, 
levitons, is analyzed in both time and energy domains. 
The conclusion is given in Sec. IV. Some details of calcu¬ 
lations are presented in Appendices A - F. 

II. EXCESS CORRELATION FUNCTION 

Let us consider an electronic one-dimensional chiral 
(or ballistic) waveguide connected to an electronic reser¬ 
voir in equilibrium, a metallic contact at Fermi energy 
/i, which serves as a source of electrons. The first-order 
correlation function for spinless electrons in a waveguide 
is defined as follows, ^^^^(1;2) = ('k^(l)^(2)), where 
^(j) = {xj,tj) is an electron field operator in the 
second quantization evaluated at point Xj and time tj, 
j = 1,2. The quantum-statistical average (...) is taken 
over the equilibrium state of the electronic reservoir. 

If a voltage V(t) is applied on the contact, then a non¬ 
equilibrium flux of particles is injected from a contact 
into a waveguide. The excess first-order correlation func¬ 
tion is defined as the difference of correlation functions 
with the voltage on and off (the subscripts “V” and “0”, 
respectively): = Qy^ — The correlation 

function G*^^^ characterizes a flux of particles injected 
from a driven contact into a waveguide on the top of the 
Fermi sea with chemical potential which plays a role of 
a fermionic vacuum. Importantly, the contribution of in¬ 
jected particles is completely separated from that of the 
underlying Fermi sea.®® 

At zero temperatures the disturbance of the vacuum 
due to the injected flux consists in the appearance of par¬ 
ticles with energy larger than /r and in the disappearance 
of some particles with energy less than In the former 
case one can speak about injection of (quasi-)electrons, 


while in the latter case about injection of holes. Below I 
will focus on an electron injection. However the presence 
of holes can be also taken into account by adding a hole 
correlation function.®®’®®’®® 

Here I am interested in the case when the amplitude 
of the voltage is relatively small, eV fi, and, therefore, 
the spectrum of injected electrons in a waveguide can 
be linearized close to the Fermi energy /r. In this case 
the excess correlation function depends on a combined 
coordinate, tj = tj — xj/v^ with the Fermi velocity, 
rather than on Xj and tj separately. 

A. A pure A-electron state 

The single-particle correlation function for a pure N- 
electron state (with possibly N oo) can be represented 
as follows, see e.g. Ref. 40, 

N 

Gw (G;t2) = (1) 

k^l 

where the mode wave functions ipk are orthonormal, 

d#fei(O'0fe2W = 4ifc2- (2) 

However, the decomposition presented in Eq. (I) is not 
always known a priori, especially if the correlation func¬ 
tion is determined experimentally. In this case the fol¬ 
lowing property of the correlation function is useful to 
prove that an electronic state in question is a pure state, 

J (3) 

As it is pointed out in Ref. 41, an oscillating poten¬ 

tial acting on non-interacting electrons of the Fermi sea 
generates a pure state. The state emitted by a single¬ 
electron source without intrinsic dephasing®®’^® is also 
expected to be pure. 

The higher-order (excess) correlation function is ex¬ 
pressed in terms of the first-order correlation function 
via the Slater determinant, 

/2, ■ • •, /mGm, *M-1, • ■ ■ , G) = 

( 4 ) 

/ G«(/i;zi) ... G(i)(/i;zm) \ 

= det 

Vg«(/m;*i) ... 

where / and i stand for the final and initial coordinates, 
respectively. The derivation of the second-order excess 
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correlation function is presented in Ref. 43. The higher- 
order correlation functions can be calculated along the 
same lines. 

Remarkably, for a fermionic A^-particle state the cor¬ 
relation functions of order higher than N all nullify, (see 
Appendix A) 


It is convenient to parametrize E by the Floquet 
(quasi-)energy, 0 < e < and by the Floquet band 

number, n, as follows E = + e + nMl = e„. With 

this notation the inverse Fourier transformations can be 
written as follows. 


f2, • ■ • ) /m; iM, *M-i) • ■ •) *i) — 0) M > N. 

(5) 

Therefore, for a single-electron state = 0, which 

is a consequence of the fact that the first-order correla¬ 
tion function is factorizable, that is the sum in Eq. (1) 
contains only one term. This property can be used to 
demonstrate a single-electron state. 


B. The energy representation 

Sometimes the energy representation is more conve¬ 
nient than the time representation. Especially if both 
electrons and holes are present. 

In general in the energy representation the two-time 
quantities like G*-^^ (ti;t 2 ) become dependent on two en¬ 
ergies. The energy and time representations are related 
via a corresponding Fourier transformation. However 
in the case of a system driven periodically with period 
To = 27r/H, the difference of two energies is not arbi¬ 
trary but can only be a multiple of the energy quantum 
hfl. It is convenient to account this property directly 
and denote one energy as E and the other one as, say, 
Em = E + mhn with m an integer. I introduce two corre¬ 
lation functions in the energy representation, and 

depending on which time <2 or G, respectively, 
can be viewed as conjugate to E. The other time can be 
viewed is conjugate to Em- The formal definitions are 
the following, 

/ OO 

-OO 

(6a) 

X - At), 

Jo Jo 


/ OO 

-OO 

(6b) 

X ^G™^‘=GW(At + t2;G), 

To 10 


where At = ti — t 2 - Note at m = 0 both functions above 
are the same, G(i)’“(E;£;) = G(^)’°“*(E; A). 


To * „=_oo 

(6c) 

OO 


G«(ti;t2) = e*^T f 

To 
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(6d) 


X ^ e-™^*^G(i)’°“‘(e„;e„+„), 


Since by its definition the correlation function satis¬ 
fies the following relation, G^^^(ti;t 2 ) = G(^^(t 2 ;G) , 

one can get for its Fourier transforms the following, 
G^^)’®"(e„; Ej,) = G^^)’°“‘(ei,; e„). The latter allows to 

represent equations entirely either in terms of gF).*" or 
in terms of G(i)-°“‘. 

The use of two different Fourier transforms is espe¬ 
cially convenient to calculate convolutions, like in the 
purity condition, Eq. (3). Applying Eourier transforma¬ 
tions presented in Eqs. (6a) and (6b) to Eq. (3) we arrive 
at the following relations specific for a pure electronic 
state. 


(7) 


GW’™(e„;e,) = GW’°"*(e„;e,), 

OO 

GW’™(e„;e,)=^;^ ^ em)&^^’^^(em; e,). 

m— — oo 

The details of calculations are given in Appendix B. 


C. Correlation function for adiabatically emitted 
electrons 

The excess correlation function can be expressed in 
terms of the Floquet scattering amplitude"*® of a peri¬ 
odically driven source, which depends on two energies. 
However in the limit of an adiabatic (slow) drive, the 
source can be characterized by the frozen scattering am¬ 
plitude, S{E,t), which depends on one energy and time. 
The dependence on energy E is the same as for the scat¬ 
tering amplitude of a stationary source. The dependence 
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on time t stems from the fact that some parameters of 
a source (and hence of a stationary scattering matrix) 
are affected by the drive. Since the drive is periodic in 
time, the scattering amplitude is also periodic in time, 
S{E,t) = S{E,t + To). At zero temperature what mat¬ 
ters is the scattering amplitude of a source evaluated for 
electrons with Fermi energy, E = ^. 1 denote such an 
amplitude as S(t) = 

In the case of a quantum-dot based source, for 
instance, a mesoscopic capacitor side-attached to a 
waveguide,the regime of a drive is adiabatic if 
the dwell time (the time an electron spends within the 
source) is small compared to the time extent of an emit¬ 
ted wave-packet.^®’^® If an electron stream is generated 
by applying an ac voltage V(t) to a metallic contact as in 
Refs. 1 and 8, then the regime is adiabatic if eV <C fi. 
The corresponding scattering amplitude is the following, 

S{t) = 50-52 

At zero temperature and in the adiabatic driving 
regime, the excess first-order correlation function is ex¬ 
pressed in terms of a scattering amplitude of the source 
S{t) as follows,®®’®® 




TtT ti — t2 


( 10 ) 


This relation, valid for a single-channel case, greatly 
simplifies an experimental set-up necessary to find a 
correlation function. Indeed, a time resolved mea¬ 
surement of the current generated by the ac source, 
= {—ie/2TT)SdS *not an interferomet¬ 
ric current,®®’®® is sufficient to access since 

G^^^t'jt) = I{t)/e. If nevertheless G^^'>{ti]t 2 ) as a func¬ 
tion of both its arguments is experimentally available, 
then Eq. (10) can be used to ensure that the driven sys¬ 
tem is in the adiabatic regime. For non-adiabatically 
driven systems the correlation function is related to 
the Floquet scattering matrix rather than to the frozen 
scattering amplitude of the source.®® As a result, equa¬ 
tion (10) is not expected to hold. 


III. A PERIODIC STREAM OF OVERLAPPING 
LEVITONS 




J_ehii-*2)tG-(i)(G;t2), 

Vn 


G'^^\ti]t2) 


l l-S*{h)S{t2) 

27r ti — t2 


( 8 ) 


The phase factor, , and the constant factor 1/v^^ 

are not important for further discussion, therefore, I will 
work below with the envelope correlation function G^®^ 
only. Note the absence of a factor in Eqs. (3), (7) 
written in terms of G*^®^ not G^®). In addition in Eqs. (6) 
written in terms of G*-®^ we have to put ^ — 0. 

The correlation function G^®^(<i;t 2 ) is not periodic in 
its arguments. However, since the scattering amplitude 
S{t) is periodic, then, as it follows from Eq. (8), the re¬ 
duced correlation function. 


9 — (G — ^ 2 ) G*-®^ (G; ^ 2 ), (9) 


is periodic in each of its arguments. Therefore, it is suf¬ 
ficient to measure G*^®^ over one period in each of its 
arguments and then to use the periodicity of g (G; < 2 ) to 
extend G^®^ to larger arguments. 

Moreover, as it follows from Eq. (8), only the diago¬ 
nal part of the correlation function contains all essential 
information. To show this let us take the limit G —>■ G 
in Eq. (8) and find G^®) {t;t) = {—i/2TT)S[t)dS*{t)/dt. 
Since the scattering amplitude is unitary {S*S = 1), 
we can represent it as S{t) = and obtain (j){t) = 

—Ttt dt'G^^^t']t'). Einally we rewrite Eq. (8) in terms 
of the correlation function only. 


In Refs. 1 and 8 the train of wave-packets carrying 
each an elementary charge was generated by applying a 
periodic sequence of the Lorentzian voltage pulses with 
integer flux, 


00 

eVstit) = 

n——oo 


2nTr 

it-nTof+n' 


( 11 ) 


to a metallic contact, as it was suggested in Refs. 20 and 
21. These elementary excitations were named levitons. 
In the equation above F,- is the half-duration of a single 
pulse and To = 2tt/Q is a period, i.e., the time delay 
between two subsequent pulses. The subscript “sf” is 
used to differentiate quantities related to a stream. For 
such a potential the corresponding scattering amplitude 
Sst{t) reads,®® 


0 ° sin ITT *+G®"" 

Sstit)= n SLit + kTo) = — 
k— — oo 

where S'L(t) = {t + iTr) / {t — iTr) is a scattering ampli¬ 
tude corresponding to a single voltage pulse, which ex¬ 
cites a single leviton. Using the equation above in Eq. (8) 
we obtain a correlation function. 



gIP (fi;G) 


* 1 - nr=-oo SI (ti + kTo) Sl (G + kTo) 
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The second line is given just for the further reference. 
However for subsequent calculations I need only the first 
line. 


A. Decomposition of the correlation function 

To bring Eq. (13) into the form of Eq. (1) let us proceed 
as follows. First, it is easy to verify the following identity, 

OO 

1 - iC/c = (1 - a^i) + (1 - X2)xi + (1 - X:i)xiX2 + ■■■ 

(14) 

OO n—1 

— ^ ^ (1 ^n) ^k' 

n—1 k—1 


Of course, the representation given in Eq. (16a) is not 
unique, since the basis single-particle wave functions 'I'n 
given in Eq. (16b) can be chosen in many different ways. 
The unitary transformations relating different bases can 
be time-dependent. 

Note that the train of hole-like levitons is excited by 
the potential of an opposite sign, Vll{t) = —Vst{t), where 
the superscript “h” denotes quantities related to holes. 
Formally it can be accounted by the following replace¬ 
ment, Ft —> —r,-, which also implies S'i(t) = [^'^(t)]* 
and 'k^(t) = ['l'i(t)]*. As a result the analogue of 
expansion given in Eq. (16a) acquires a minus sign, 

(ti;t 2 ) = -Er=-oo [^n(^i)]*^nfe)- The minus 
sign is a mere consequence of the fact that we consider 
an excess correlation function. 


B. Periodicity properties 


Now one can see that the next identity is also correct. 


OO OO OO 

1 - Xfe = (1 - 2 : 1 ) Xfc -f (1 - X2) Xk 

k^l k^2 k^3 

(15) 

OO 0000 

+ (1 - 2 : 3 ) n H- = '^{^-Xn) n 

fc—4 n—1 fc—n+1 

The equation above is also valid if the lower limit is minus 
infinite, k —> — 00 . 

Then let us use Eq. (15) with an infinite lower limit 
and Xk = S^{ti + kTo)SL{t 2 + kTo), and rewrite Eq. (13) 
(its first line) as follows. 


OO 

{ti;t 2 ) = ^n(^i)^n(i 2 ), (16a) 

n——(x> 


where the single-particle (envelope) wave functions. 


OO 

T„(t) = 4'i,(t-I-nTB) Skit + kTo), (16b) 

k—n-\-l 


are orthonormal, 





= Sr, 


(16c) 


Here ^'^(t) = v^TV TT/ {t — zTt-) is an envelope wave 
function of a single leviton in a one-dimensional 
channel.Note that the full leviton’s wave function is 
4'i,(<). To get Eq. (16a) we also 
need to use the following relation, 1 —S’l{ti+nTo)SL{t 2 + 
nTo) = - 27 ri(ti - t2)'^l{ti + nTo)'^L{t2 + nTo). 


The basis wave functions 4'„(t) in Eq. (16b) are not 
periodic in time. Therefore, the correlation function 
G[Wti;t 2 ), Eq. (16a), is also not periodic in its ar¬ 
guments. However, the reduced correlation function, 
gstiti;t 2 ) = [ti -t 2 )Gii^ {ti;t 2 ), Eq. (9), is periodic in 
each of its arguments, see Fig. 2. 

The wave functions possess the following discrete 
time-translation symmetry. 


«'„(< + ro) = ^„+i(t), (17) 


which reflects the periodicity of a particle flux. 





FIG. 2. (Color online) The real part of the reduced first-order 
correlation function of the stream of levitons, gst{ti',t 2 ) = 
(h — t2) Gj)' (ti; 12 ), is given as a function of its arguments t\ 
and t 2 normalized to the period To- For G^)^ see Eq. (16a). 
The width of a leviton is chosen to be 2TrlTo = 0.1. 
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The correlation function at coincident times is peri¬ 
odic, 


GW(t;t)=G«(t + r;t + r)= ^ .(18) 

n— — oo 

Using Eq. (17) one can then evaluate the mean number 
of electrons (levitons) per period, 


To 

{N,t)= f ”rftG«(t;i) 

Jo 


rTo 


dt 


OO 




/ OO 

dt\^L{t)\^ = 1, 

-OO 


(19) 


Finally using Eq. (10) the full correlation function of a 
single leviton, G^^^ (ti;t 2 ), can be reconstructed. 


C. The purity condition in time and energy 
representations 

From Eqs. (16) it immediately follows that Eq. (3) 
(without a factor u^) holds, hence the state is pure. 
Therefore, the levitons are created in a deterministic 
fashion. 

The purity of a state can also be inferred from a corre¬ 
lation function in the energy representation. Substitut¬ 
ing Eq. (16a) with wave functions from Eq. (16b) into 
Eqs. (6b) and (6a) we calculate. 


which is one, independently of the ratio Ft-/T o, which 
characterises the degree of overlap between the levitons. 

Note that the periodicity of a flux allows to express 
the correlation function of a single leviton, G^^^ (U; ^ 2 ) = 
d'^(ti)'I'i(t 2 ), in terms of an experimentally accessi¬ 
ble correlation function of a stream, GgJ^(<i;< 2 ). For 
this purpose we use Eqs. (17), (18) and relate a dis¬ 
crete Fourier transformation, of a periodic func¬ 

tion Ggj^ (t; t) to a continuous Fourier transformation, 
of a non-periodic function G^^^ (t;t), as follows. 


Gil\k)= 

Jo /o 


0^0 


2 


= E 


n= —00 

0° -(n-l-l)To j-Ll 


n——oo 
roo j.f 


nTo 


To 


r %Gf = ^Gf{u: = kn), 

J-oo >0 /o 

( 20 ) 


GW’“(e„ > 0; e, > 0) = GiJ)’°“‘(e„ > 0; e, > 0) 

( 22 ) 

= 2sinh(Ur^)e-^^-(”+'^+F. 


Remember that e„ = e -I- nhVL with the Floquet energy 
0 < e < Ml. Therefore, e„ > 0 implies n > 0. For other 
energies the correlation function is zero, what is specihc 
for purely electronic excitations.^® Note that G^^^ (k) = 
(l/ro)exp(-|A:|Ur,). ^ 

By applying the inverse Fourier transformation, 
Eq. (6d) or (6c), to Eq. (22) we arrive at G^^^ {ti\t 2 ) 
in the form of the second line of Eq. (13), as it should be 
(see Appendix C). The two functions above, and 

g(i), out, satisfy Eq. (7) (without a factor v^) charac¬ 
teristic for pure states. Note that in the present case the 
sum over m in Eq. (7) runs from m = 0. 

Thus, the purity of a state of the stream of levitons is 
demonstrated in both time and energy representations. 


D. Factorization in the time representation 


where I made a substitution t' = t — nTo and took 
into account that Tokl = 27r, hence exp(—ifcUnTo) = 1- 
Note that a single-time Fourier transformation over t = 
h = M G(i)(fc), is related to a two-time Fourier trans¬ 
formation introduced in Eq. (6b) as follows, G^^\k) = 
(l/ro)E“-ooG«’“ (enUn+fe). 

Equation (20) defines G^l\uj) only on a discrete set of 
frequencies, w = kkl. However, in the regime of adiabatic 
emission the Fourier coefficients G^^^ (w) change smoothly 
on the scale of H. This allows us to find g''^\ui) at other 
frequencies by interpolation and calculate. 



The sum in Eq. (16a) involves more than one term, 
therefore, in general the state of the stream of levitons 
is a multi-particle state. As a result, all higher-order 
correlations functions are non zero. The multi-particle 
state is formed due to an overlap of the levitons excited 
during different periods, which is the case when Ft- ~ Tq. 

In contrast, if the levitons do not overlap, which is the 
case when Ft- <C To , then the state of a flux is the direct 
product of states of individual levitons. This formally 
follows from Eq. (16b), where limr^/ 7 j_s.o SL{t + kTo) = 1 
and, therefore, 'I'„(t) = d>L(t -|- nTo). The state with a 
wave function exists only at a period for which 

t « —nTo and it vanishes at other periods. 

In practice whether the states of electrons in the stream 
are single-particle states or not can be verified by ex¬ 
amining whether the correlation function, G*-^^ (ti;t 2 ) is 
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factorized into two factors each of which depends only 
on one argument or not. If the first-order correlation 
function is factorizable then, as I already mentioned, all 
higher order correlation functions are zero, in particular, 
= 0. Therefore, with increase of the ratio 7o/(2rT-) 
the second order correlation function is expected to de¬ 
crease. To confirm this expectation let us evaluate the 
mean number of pairs of levitons per period, (IV 2 ), 

{N2) = (23) 


To/i^Tr) —t 00 , see Fig. 3 (a blue dashed line). The 
details of calculations of Pi are given in Appendix D. 

E. Mean energy per leviton 

Another indication that the overlapping levitons do af¬ 
fect each other is an increase of the mean energy per 
leviton, (e), with increasing overlap. Using a single¬ 
particle distribution function for the stream of levitons, 
fstie) = Git^’“(e„;e„) = /st(n), Eq. (22), one can calcu¬ 
late, 


which should decrease with decreasing the second-order 
correlation function G),f . Representing G).^ in terms 
of according to Eq. (4) and taking into account 
Eq. (19), one can express {N 2 ) in terms of an experi¬ 
mentally accessibly first-order correlation function. 


{N 2 ) = 1 - 



dtdt' 


^( 1 ) 




(24) 


where G^J^ is given in Eq. (13) [or in Eq. (16a)]. The 
dependence of {N 2 ) on the ratio 7o/(2r,-) is shown in 
Eig. 3 (a red solid line). The number of pairs of levitons, 
which can be detected during a period, gradually de¬ 
creases with increasing the distance between the levitons, 
which demonstrates an evolution from a multi-particle 
state to a single-particle state. At the same time the 
probability to measure exactly one leviton during a pe¬ 
riod, Pi, gradually increases up to one in the limit of 



FIG. 3. (Color online) The mean number of pairs of levitons 
detected during a period, {N 2 ), Eq. (24), (a red solid line) 
and the probability to measure exactly one leviton during 
a period. Pi, Eq. (D2), (a blue dashed line) are given as a 
function of the ratio of the time delay between two subsequent 
levitons, To, and the width of a leviton, 2rT. Notice the 
different axes used for {N 2 } and Pi, which are indicated by 
horizontal arrows. In the limit when levitons do not overlap, 
7o/{2Tt) S> 1, we have {N 2 } — >■ 0 and Pi — >■ 1, which means 
that the multi-particle state is not formed. 


(e) 


/ dee fSt (e) 
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hQ.^{n + Q.h)fst{n) 

n—Q 


hn 


coth (nr,-) 
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2r, ’ To ^ ^ 
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(25) 


The energy (e) gradually increases with increasing over¬ 
lap from the value characteristic for an iso¬ 

lated leviton to the value hn/2 characteristic for a dc 
bias The latter is what Ust(f), Eq. (11), looks like at 
r.» To.^^ 


F. Accidental factorization in the energy 
representations 

In contrast to a correlation function in the time repre¬ 
sentation, the correlation function of the stream of levi¬ 
tons in the energy representation, Eq. (22), is factorized 
at any degree of overlap, Ggj^’“(e„; Cj.) = ipst(n)ipst{v), 
where (pst{n) = •\/2 sinh (nFr) . Such a fac¬ 

torization results in vanishing of a two-particle distribu¬ 
tion function (see Appendix E), fst{^n,^v) — 0. This 
result can be interpreted as a zero probability to mea¬ 
sure two levitons with fixed energies, e„ and However 
this fact does not deny a possibility to measure two levi¬ 
tons during a time period of duration To, see Fig. 3. The 
point is that a leviton has no a definite energy and its 
projection into the state with a definite energy can only 
be done via measurement over an infinite time. The pos¬ 
sibility to measure two levitons during a finite time is 
also clearly demonstrated by the theory of waiting time 
distributions.®® 

The factorization mentioned above is to some extend 
accidental and specific to a driving protocol given in 
Eq. (11). For a general driving protocol the correlation 
function in the energy representation reveals the presence 
of multi-particle correlations. The energy representation 
is particularly useful to demonstrate the presence of cor¬ 
relations between electrons and holes created by an ac 
voltage. As an example the case of a harmonic driving 
voltage is shortly discussed in Appendix F. 










IV. CONCLUSION 


The excess first-order correlation function allows us to 
differentiate unambiguously whether the state of a peri¬ 
odic stream of electrons injected into a waveguide is a 
multi-particle state or can be represented as the product 
of single-particle states existing during each period sepa¬ 
rately. In the latter case the second- and all higher-order 
correlation functions, which all are expressed in terms of 
the first-order correlation function, vanish while in the 
former case they do not. 

Motivated by recent experiments^’® I analysed a corre¬ 
lation function of a train of single electrons, levitons, ex¬ 
cited by a periodic sequence of Lorentzian voltage pulses 
with a unit flux each. I decomposed this function into the 
sum of single-electron contributions. This allows us, first, 
to verify that the state of a stream is pure rather than 
mixed. Second, such a decomposition in the time rep¬ 
resentation allows us to show clearly that decreasing the 
width of a leviton makes the state of the electronic stream 
evolving from a multi-particle state to a single-particle 


state. In the latter case the correlation function is fac- 
torizable while in the former case it is not. Unexpectedly 
I found that in the energy representation the correlation 
function is factorizable at any width of a leviton. This 
hinders a direct use of energy-resolved measurements to 
differentiate single- and multi-particle states in the case 
of a stream of levitons. The possible reason is that such 
measurements involve averaging over an infinite time in¬ 
terval, over which the single-particle wave functions are 
orthogonal due to the Fermi statistics irrespective of the 
degree of overlap.Due to orthogonality even the over¬ 
lapping particles contribute independently. 
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Appendix A: {N -|- ... )th order correlation function for an V-electron state 

The determinant is calculated as follows. 


/ An . 

• Aim \ 

M 

M 





■ -^MjM ~ 'y ^ Ajil . . . Ajjj^m- 

(Al) 

^ Ami • 

• Amm / 

ii>'" om=i 

jl.'" .JM = 1 



Here is the Levi-Civita tensor with the following properties: (i) ei 2 ... 7 v = I, (h) if any two indices are the 

same, then = 0; and, finally, (iii) when the two any indices are interchanged then the tensor changes a sign, 

e...ii...i 2 ... = —e...i 2 ...ii...- Using Eq. (Al) in Eq. (4) with M ^ N we get, 

M N 

As it follows from Eq. (I), there are N different functions '0^, k = 1,2,..., N. Since M > N, then for any term in 
equation above among M factors (- = 1,2,..., M, one can always find at least two factors, say for € = a and 

£ = /3, for which = kp. Let us take a close look at the product of such two factors (we also keep the Levi-Civita 
tensor). 


= -e...j,3...jc...V’L(/“)^fcc{*u)^L(//3)V'fec(Uo)- (A3) 

Here we utilized the fact that ka = kp but in general /„ 0 fp and ^ 0-^. Going from the second to the 
third line we made an interchange ja ^ jp and took into account that the Levy-Civita tensor is antisymmetric, 
- The interchange ja ^ jp does not change Eq. (A2), which involves sums over all Jq’s. On 
the other hand, according to Eq. (A3), such an interchange should reverse a sign. These two consequences of the 
interchange ja ^ jp are consistent if and only if the entire sum is zero. Therefore, Eq. (5) is proven. 
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Appendix B: The derivation of the purity condition in the energy representation, Eq. (7) 

Let us rewrite Eq. (3) in terms of the Fourier coefficients of a correlation function. It is convenient to transform 
according to Eq. (6d) and according to Eq. (6c). As a result we arrive at the following, 


Vfj, j dt f 
J — oo J 0 


t -t °° °° 

) gmO(ti-t) "y g-zmPt 


(Bl) 


n——oo 


m— — oo 


rh^ 7 / oo 

de t-t 2 


= ^ V. 


^ 

n' = —OO m' — — oo 

l^hQ. 1 oo oo oo 

^ n— — oo m—— oo m'——oo 


where At = ti—^ 2 - Note that the integral over t produces h5 (—e — (n + m)Ml + e' + {n' + Since both Floquet 

energies e' and e belong to the same interval of length this delta function implies e' = e and n + m = n' + m'. In 
the last line of Eq. (Bl) we additionally changed n + m ^ n, m ^ —m, and m' —^ —m!. 

To proceed let us perform the Fourier transformation on the left hand side either for in- or for out- representation 
using Eqs. (6a) or (6b), respectively. First we use Eq. (6a), put t 2 = ti — At, and find 


/ oo pTo 

-(50 J 0 


dt 

0 "To 




1 / oo oo oo 


'0 


n— — oo m— — oo m——oo 




(B2) 


m — — OO 


where i, v are integers. Next we use Eq. (6b), put ti = At -|- t 2 , and arrive at the following, 

rl+ phQ I f ^ 

-oo Jo 10 Jo h ^ 
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(B3) 


m— — oo 


So from Eqs. (B2) and (B3) we get, 

OO 

G(i)’“(e„;e,)=n^ ^ G«’°“*(e„; e™)G«’“(e™; e,), 


G«’°“‘(e„;e.)=t^M E G(i)’°“*(e„;e™)G(i)’“(e™;e,). 


(B4) 


m— — oo 


Since the right hand sides of these equations are the same, the left hand sides are also the same. Therefore, we arrive 
at Eq. (7). 


Appendix C: The inverse Fourier transformation of the correlation function of the stream of levitons 


Let us apply the inverse Fourier transformation, say, Eq. (6d), to Eq. (22), where only n > 0 and v = m + n > 0 
do contribute. As a result we obtain the following equation, 


hQ 7^ oo 

de 


GiP(G;t2)=2sinh(L!r.) / 
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(Cl) 
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which coincides with the second line of Eq. (13). 
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Appendix D: Probability to measure exactly one leviton per period 


The probability, to detect a leviton with a wave function Eq. (16b), during the time interval (0;7o) is the 
following. 


rTa 


Pn= dt|4'„(t)| = 


r*(n+l)7o 


/nTo 


(t; t) = 


arctan - arctan 


(Dl) 


The sum of all p„’s gives the mean number of levitons detected during a period, = {^st) = 1, which 

agrees with Eq. (19). However the actual number of detected levitons varies from period to period. In particular, the 
probability to measure exactly one leviton per period irrespective of its state is given by the following equation, 

OO OO OO 

^1= Y. p- n (i-Pfc)=Po E (d 2 ) 

n—— OO k^n,k——oo n— — oo P^ 


where Pq = n^-oo(l ~ Pk) is a probability to do not measure any leviton. In the equation above the probability pn 
to detect a leviton with a wave function is weighted by the product of probabilities, {1 — pk), to do not detect a 
leviton in any other state, k ^ n. Note that Pi = {N^t) = 1 if and only if all pk = 0 but some = 1. This is the 
case in the absence of overlaps, 7o/(2r,-) —)■ oo, which is clearly seen from Fig. 3 (a blue dashed line). In this limit 
the number of particles detected during a given period does not fluctuate. 


Appendix E: Two-particle distribution function 

The two-particle distribution function is defined as follows, 

rlk fT) roo /■oo , 

f{en,en+m)= / 7 ^ 7 ^ / dr'e“*''"+"‘^G'(^)(T-t t,r'-f t';<',t) 

Jo '0 Jo >0 J— oo J— oo 

— f (^n)/ (^n+m) H“ /(^n? ^n+m) ; 


(El) 


where /(e„) = (e„; e„) = (e„; e„) is a single-particle distribution function. The irreducible part reads, 

/■T) Pk f'P’ rtf' r°° r°° ' 

<5^/(en,e„+r„) = (-1) / 7=^ 7=r j t;t')G(^)(r'-f t';t). (E2) 

Jo ^0 Jo ^0 J— oo J— oo 


After the following shifts, t ^ t — t + t' and t' ^ t' — t' +1, the equation above becomes 
<5^/(en,en-Hm) = (-1) / ^6 


To pp , r°° fPo Pk / 

—OO 


/o 'To 

= (_l)G(i)'°“‘(e„; e„+^)G(i)>°“‘(e„+,„; e„), 


/o '0 


To 


(E3) 


where we used Eq. (6b) written in terms of G*^^^ not G*-^^. Substituting the equation above into Eqs. (El) we arrive 
at the following equation. 


/(e„, e,) = G(i)’°“*(e„; e„)G(i)’°“*(e,; e,) - G(i)’°“‘(e„; e,)G(i)-™‘(e,; ej, 


(E4) 


which is used in Sec. IIIF. 

Appendix F: First-order correlation function of an 
electron-hole flux generated by an AC voltage 

In contrast to a periodic sequence of Lorentzian volt¬ 
age pulses, Eq. (II), the general ac voltage excites both 


electrons and holes. As an example let us consider a 
harmonic voltage, I4c(t) = E cos (fit), applied across a 
single-channel ballistic electronic waveguide connected to 
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two metallic contacts. The corresponding scattering am¬ 
plitude reads (up to an unimportant constant phase fac¬ 
tor), Sac{t) = exp {—*1^ sin(nt)}. Without the loss of 
generality one can put et^ > 0 and use the following ex¬ 
pansion, 


n——oo ^ ' 

where Jn is the Bessel function of the first kind and of 
the integer order n. Then we represent the correlation 
function, see Eq. (8), as follows. 
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where we also used the following identity, 


(F2) 


^ ' Jn+k ( '^’n.+m \ “ ^k,m- (F3) 

n— — oo 
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Then we split the sum over n in Eq. (F2) into the two 
parts, for n > 0 and n < 0, and arrive at the following. 


GW(ti;i2)= E E E J> 

m— — oo n=0 k—n-\-l 

= Gi‘> (A() Jt ^ ^ e 


with 


r{n+ 

Gi') (At) = / 

J nhQ 


(n+l)/in 


(F5) 


The correlation functions Gn'^ for different n are mu¬ 
tually orthogonal and do satisfy the purity condition, 
Eq. (3). Therefore, the corresponding states constitute a 
convenient basis for representing the state generated by 
a harmonic voltage.®® 

As I already mentioned, electron and hole contribu¬ 
tions to the correlation function have different signs, 
which hinders a direct use of Eq. (3) for the state purity 
verification. To assess electronic and hole contributions 
separately it is useful to go over from a time domain to 
an energy domain. After performing the Fourier trans¬ 
formation according to Eqs. (6a), (6b) we have, 


-f — QW,out/ 'I — 

'-^ac — '-'ac — 


' OO 

Jk+n (l^) Jk+u (l^) , ^ > 0, 

— oo 

- Z) (l^) Jk+u (l^) , n < -1. 

„ k—0 


(F6) 


Here both n and v are integers. Remember that e„ = 
e + nhO. and the Floquet energy e is chosen to be positive, 
0 < e < Ml. 

If in the equation above both energies belong to ei¬ 
ther electron, e„,ei, > 0, or hole, e„,ey < 0, sectors the 
equation (7) holds, what tells us that the emitted state 
is pure. Note that the sum on the right hand side (RHS) 
of the second line of Eq. (7) comprises both purely elec¬ 
tron (hole) and electron-hole contributions. Therefore, 
both an electron coherence and an electron-hole coher¬ 
ence are equally important to maintain the entire state 
pure. A non-zero value of CiZ*” in the electron-hole 
sector, e„ > 0,ey <0 or e„ < > 0, indicates the 

presence of electron-hole correlations.®^’®®’®® Note that 
in the electron-hole sector the sum on the RHS of the 
second line of Eq. (7) is zero due to cancellation of elec¬ 
tron and hole contributions, which is specific for states 
with equal number of electrons and holes. 
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